METHOD FOR MEASURING GRADIENT INDEX DISTRIBUTION OF ROD LENS 



BACKGROUND OF THE INVENTION 

The present invention relates to a method for measuring 
a gradient index distribution of a rod lens and particularly 
to a method for measuring a gradient index distribution of a 
gradient index rod lens by calculating higher-order index 
distribution coefficients of the rod lens on the basis of 
measurement of curvature of field. The method according to the 
invention is a technique particularly useful for evaluation of 
optical performance of a small-diameter rod lens. 

As known commonly, a gradient index rod lens is a lens having 
a columnar transparent body to which a distribution of refractive 
index symmetric with respect to the optical axis of the lens 
is given. The refractive index of the lens is distributed so 
that the refractive index is high on the optical axis but is 
reduced continuously toward the periphery of the lens. The 
gradient index rod lens of this type has been used as a collimator 
lens or the like in an optical communication system, an optical 
measurement control system or the like because reduction in size 
and weight can be achieved. In addition, a lens array constituted 
by a large number of the gradient index rod lenses of this type 
arranged regularly in the form of an array has been used as a 
scanning optical system in a copying machine, a facsimile machine, 



a printer or the like. 



Various methods have been proposed for giving a gradient 
index distribution to a columnar transparent body (glass rod) . 
The method put into most practical use is an ion exchange method. 
This is a method in which a glass rod containing 
high-refractive-index ions is immersed in molten salt containing 
low-refractive-index ions to disperse the two types of ions into 
each other to thereby form a distribution of refractive index 
(gradient index distribution) approximately proportional to the 
distribution of ion concentrations. 

The optical performance of the rod lens of this type mainly 
depends on the shape of the gradient index distribution, it 
is therefore necessary to control the distribution for production 
of the lens. For this reason, the gradient index distribution 
is required to be measured accurately, m addition, index 
distribution coefficients are very important as basic data for 
evaluation of variation in ion exchange and lens design or for 
system design using such a rod lens. 

As a method for obtaining a gradient index distribution 
of a gradient index rod lens, there has been heretofore used 
a method of calculating back the gradient index distribution 
by measuring spherical aberration of a P/4 lens (in which P 
representsaparaxialperiodlength) (see "Measurement andAnalysis 
of Aberration of Gradient Index Lens", Optics Vol.11, No . 6 
(December 1982) ) . 



t 



Also in the rod lens, spherical aberration can be obtained 
by direct measurement of the locus of laser light. That is, 
because laser light incident on one end surface of a lens to 
be inspected passes through the lens and exits from the other 
5 end surface of the lens, the locus of light rays is obtained 
by observationof exit light rays . When themeasurement is repeated 
while the position of incidence of light is changed, flux of 
exit light rays is obtained so that spherical aberration can 
be obtained. The related-art method is a method in which a light 
QO ray equation is solved by a perturbation method in consideration 
W\ of a gradient index distribution of a gradient index rod lens 
m up to higher-order terms to thereby obtain approximate solutions 
m to parallel light incidence so that spherical aberration of the 
L lens is measured by application of the approximate solutions 
LM.5 to thereby obtain index distribution constants. 

E F 

5* In the related-art method using the measurement of spherical 

z- r h 

aberration, it is however necessary that the measurement is 
repeated while the position of incidence of light is changed 
in a direction of the radius of the rod lens. Accordingly, it 

20 is difficult to measure spherical aberration particularly when 
the small-diameter rod lens has a diameter not larger than about 
1 mm. With the advance of reduction in size of various kinds 
of optical devices in recent years, the rod lens to be incorporated 
in each of the optical devices is required to have a further 

25 smaller diameter. Accordingly, the index distribution 
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coefficients can be hardly obtained by the related-art method. 

In addition, the related-art method needs laser light as 
a light source . Hence, there is also a problem that the wavelength 
for measurement is limited to the wavelength of laser light used. 
5 SUMMARY OF THE INVENTION 

An object of the invention is to provide amethod for measuring 
a gradient index distribution of a rod lens in which higher-order 
index distribution coefficients can be obtained accurately even 
if the rod lens has a small diameter and in which the wavelength 
qLO for measurement can be selected relatively freely, 
m According to the invention/ there is provided a method of 

Hi measuring a radial gradient index distribution n(r) of a rod 
S{ lens by calculating higher-order index distribution coefficients 
U indicating the gradient index distribution n(r) when n(r) is 
2l5 given by the expression: 

§ n(r) 2 = n 0 2 -{l - (g-r) 2 + h 4 (g-r) 4 + h 6 (g-r) 6 

- % 

R5 + h 8 ( g • r ) 8 + • • • } 

in which r is a radial distance measured from an optical axis, 
n 0 is a refractive index on the optical axis, % is a secondary 
20 index distribution coefficient, and h 4 , h 6 and h 8 are higher-order 
index distribution coefficients, the method comprising the steps 
of: 

(1) processing the rod lens so that the optic-axial length 
of the rod lens is approximately equal to P/2 (in which P represents 
25 a paraxial period length (pitch) defined as P = 2n/g) or 
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approximately equal to an integer multiple of P/2 and so that 
opposite end surfaces of the rod lens are shaped like parallel 
planes; 

(2) setting a patterned surface as an object surface in 
the proximity of one end surface of the rod lens and forming 
an image surface in the proximity of the other end surface of 
the rod lens by irradiating the patterned surface with condensed 
monochromatic light ; 

(3) obtaining the position of a paraxial focal point and 
the curve of curvature of field by observing the image surface; 
and 

(4) calculating back higher-order index distribution 
coefficients by a fitting process on the basis of the position 
of the paraxial focal point and the curve of curvature of field. 

Here, preferably, the patterned surface has a striped 
structure in which a large number of straight lines are arranged 
in parallel with one another or a lattice structure in which 
a large number of straight lines are arranged crosswise and in 
parallel with one another so that the positions of focal points 
on a plurality of lines are measured in a direction of the arrangement 
of lines from the center of the lens (optical axis) to thereby 
obtain the curves of curvature of field of meridional image surface 
in accordance with distances from the optical axis of the lens. 

The pre sent disclosure relates to the subject matter contained 
in Japanese patent application No . 2001-033899 (filed on February 



9, 2001), which is expressly incorporated herein by reference 

in its entirety. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is an explanatory view showing the relation between 
5 optical path and curve of curvature of field in accordance with 
the change of the value of h 4 . 

Fig . 2 is an explanatory view showing an example of a measuring 
apparatus . 

Figs. 3A and 3B are explanatory views showing examples of 

5 ;. 

Q0 the patterned surface. 

gj Fig. 4 is a flow chart showing a method for determining 

rg index distribution coefficients h 4/ h 6 and h 8 . 

a jfSj 

m Fig. 5 is an explanatory view showing fitting of the index 

n distribution coefficients. 

jjj5 Fig. 6 is an explanatory view showing another example of 

y the measuring apparatus. 

HI DESCRIPTION OF THE PREFERRED EMBODIMENT 

In the invention, typically, curves of (particularly 
meriodinal) curvature of field of a lens having a lens length 
20 of about P/2 (in which P is a pitch indicating the length of 
a paraxial snaking period in the lens) is measured so that index 
distribution coefficients of the lens are calculated back 
softwarewise to be fitted to the curve. 

As shown in Fig. 1, when an image is formed between opposite 
25 end surfaces of a P/2 rod lens, an object surface and an image 
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surface are located to be symmetric with respect to the rod lens. 

Accordingly, there is neither comatic aberration nor distortion 
caused by asymmetry, so that the image has no aberration but 
spherical aberration and (sagital and meridional) curvature of 
field. Hence, the curvature of field can be measured accurately 
if the spherical aberration of the lens is not extremely large. 

Incidentally, if a rod lens having a lens length n times 
as large as P/2 (in which n is an integer, that is, the lens 
length is IP, 1.5P, 2P, ...) is used, the quantity of aberration 
is also increased by n times so that measuring accuracy can be 
improved. It is, however, neces-sary to select the value of n 
suitably in accordance with the lens because the image becomes 
too blurred to make the measurement accurately if the lens is 
large in spherical aberration or is striated or asymmetrical. 

If a rod lens having a lens length of P/4 is used, a pattern 
is set in a distant place to make it difficult to measure a meridional 
image surface because of the influence of comatic aberration 
on the focussed image. Accordingly, the use of a lens having 
such a lens length is unsuitable. 

As an example for indicating the relation between index 
distribution coefficient and curvature of field, there was 
calculated the curvature of field of a P/2 rod lens having the 
following specifications . 

(Specifications of Model Lens) 

Effective Radius r 0 = 0.125 mm 



Optic-axial Refractive Index n 0 = 1.682 

2 value = 3.15 /mm 

Lens Length Z = 0.997 mm (P/2) 

Object Surface and Image Surface: opposite end surfaces 
of the rod lens 

Object Height: 0.125 mm 
(Optical Path Chart and Field Curvature) 

Fig. 1 shows optical path charts and curves of curvature 
of field in the respective cases of index distribution coefficient 
h 4 equal to -1.0, +0.67 and +2.0. 

In the case of h 4 = -1.0, a meridional image surface (M) 
and a sagital image surface (S) are both minus. 

In the case of h 4 = +0.67, the meridional image surface 
(M) is substantially flat but the sagital image surface (S) is 
minus . 

In the case of hu = +2.0, the meridional image surface (M) 
and the sagital image surface (S) are both plus. 

As is obvious from Fig. 1, the measurement of field curvature 
data on the basis of the meridional image surface (M) is more 
accurate because the quantity of change of the meridional image 
surface (M) relative to h 4 is larger than that of the sagital 
image surface (S) . In addition, when a striped pattern is used, 
the measurement can be made more easily on the basis of the meridional 
image surface because the meridional image surface is easy to 
view. 



Fig. 2 shows an example of a measuring apparatus. This 
apparatus mainly has a stage 10 movable finely vertically, a 
linear gauge 12 for measuring the height of the stage 10, a light 
source 14 located below the stage, and a microscope 16 located 
above the stage. A transparent parallel-striped pattern 18 is 
placed on the stage 10. A lens 20 to be inspected is placed 
on the pattern 18. The pitch of the parallel-striped pattern 
18 is selected to be about 1/20 as large as the diameter of the 
lens. The parallel-striped pattern 18 is placed on the stage 
10 so as to be striped longitudinally. The lens 20 to be inspected 
is a rod lens having a lens length of about P/2. Opposite end 
surfaces of the lens 20 to be inspected are shaped like planes 
perpendicular to the optical axis. Illumination light from the 
light source 14 is changed into monochromatic light with a 
measurement wavelength by an interference filter 22 and condensed 
onto the lower end surface of the lens 20 by a condenser lens 
24. A stop 26 is provided on the exit side of the condenser 
lens 24 so that the numerical aperture (NA) of the illumination 
light is selected to be about 0.1. The magnification of the 
microscope 16 is selected so that the whole upper end surface 
of the lens 20 to be inspected can be viewed. An image surface 
near the upper end surface of the lens 20 to be inspected is 
observed with the microscope 16. 

As shown in Fig . 6, a CCD 102 may be connected to the microscope 
16 to obtain digital data of the image obtained through the 



microscope 16, and a computing unit 100, such as a PC, may be 
connected to the CCD 102 for image processing, and for measurement 
and computation of values required for determining the higher-order 
index distribution coefficients h 4/ hr, , ha ... to be described 
later* In this case, the computing unit 100 has a storage medium 
which stores therein or which can install therein a program for 
executing a method of the present invention for determining the 
higher-order index distribution coefficients h 4 , h 5 , h 9 ... to 
be described later, and also has an input section by which 
predetermined data, known values, measured values, etc, can be 
inputted. The computing unit 100 may further be connected to 
drivers for the microscope 16, the stage 10, the light source 
14, etc. for control of these components, and also to the linear 
gage 12 for obtaining the height data of the stage 10. In Fig. 
6, reference numeral 104 designates a driver for controlling 
the microscope 16 . The computing unit 100 may further be connected 
to a driver 106 for moving the movable stage 10 in X-, Y- and 
Z directions. The computing unit 100 can adjust the height of 
the stage 10 (i.e. the position of the stage in the Z-direction) , 
and obtain focal points based on image data obtained by the CCD 
102 andheight data obtainedby the linear gage 12 . Suchmeasurement 
can be carried out at plural locations under the control of the 
computing unit 100 by moving the stage in lateral directions 
(i.e. in the X-direction and in the Y-direction) through the 
driver 106. The computing unit 100 can calculate and display 



the higher-order index distribution coefficients h A , h 6 , hg 
based on the data thus obtained. 

As for the measurement, the lens length is first measured 
accurately and the center refractive index no of the lens is 
5 measured in advance. If the magnification of a microscope 

objective lens is set to be high (that is, NA is set to be large) 
in order to measure a small-diameter lens, measuring accuracy 
may be deteriorated because spherical aberration of the lens 
is picked up so that an image is blurred or a focal point is 
30 moved. To prevent the deterioration of measuring accuracy, the 
yl NA of illumination light may be preferably set to be small, for 
g| example, about 0.1, and it is undesirable to use any diffusing 
fi\ plate or the like. If the NA of the microscope objective lens 

5 

p* is not larger than about 0.1, the consideration described above 

jif.5 is not required because light rays with large NA are cut off. 
i? Approximate values of index distribution coefficients can 

S M 
T* r r 

jFH S 

iU be calculated back by use of optical design software available 
on the market. For example, "Oslo Six" made by Sinclair Optics, 
Inc. in the U.S. may be used. 
20 Although the example has shown the case where the patterned 

surface located on the object surface has a striped pattern in 
which a large number of straight lines are arranged in parallel 
with one another at regular intervals as shown in Fig. 3A, the 
invention may be applied also to the case where the patterned 
25 surface has a lattice pattern in which a large number of straight 
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lines are arranged crosswise and in parallel with one another 
at regular intervals as shown in Fig. 3B . 

Fig. 4 shows steps for determining higher-order index 
distribution coefficients h 4 , h 6 and h 8 . 
(Preparation for Measurement) 

The lens length Z (approximately equal to P/2) of the lens 
(having opposite ends shaped like planes) to be inspected is 
measured accurately. 

The value of center refractive index n 0 of the lens is measured 
in advance. 

The lens to be inspected is placed on the patterned surface 
of the measuring apparatus shown in Fig. 2. 

The magnification of the microscope is selected so that 
the whole end surface of the lens can be viewed. 
(Measurement of Field Curvature) 

The position of a paraxial focal point Afo is measured in 
the center of the lens with the lens end surface as the origin 
of the linear gauge (in which the outside of the lens is regarded 
as being plus and the inside of the lens as being minus) . 

The positions of (meridional) focal points Afi, Af 2 , Af 3 ,.~ 
of the respective lines in accordance with distances r x , r 2 , r 3 , 
... from the optical axis are measured in the same manner as described 
above . 

(Setting of Initial Parameter with Design Software) 

The lens radius is measured on the lens to be inspected. 



The optic-axial refractive index n 0 is measured on the lens 
to be inspected. 

The lens length Z is measured on the lens to be inspected. 

The index distribution coefficient cj value is initialized 
to n/Z. 

h 4 , hs and h 8 are all set to be zero. 

The object surface is made coincident with one end surface 
of the lens. 

The image surface is set in a position far by the paraxial 
focal point position Af 0 from the other end surface of the lens. 

The light source with respect -to the obj ect surface is provided 
as a telecentric structure (in which main light rays are parallel 
with the optical axis) . 
(Accurate Determination of c[ Value) 

The cj value is adjusted so finely that the paraxial focal 
point of the lens to be inspected is made coincident with the 
image surface. 

(Setting of Object Height and Meridional Focal Point Position) 
Points of object height r x , r 2 , r 3f ... are set on the object 
surface . 

A distance Zi from the image surface = Afi - Af 0 is set as 
a target value of the meridional focal point position of luminous 
flux exiting from the point of object height r x . 

z 2 , z 3 , ... are set in the same manner as described above. 
(Definition of Optimizing Function) 



The difference from the target value is defined as follows: 
Azi = zi' - Zi 

in which z x ' , ZJ , Z 3 ' , ... are calculated values of the meridional 
focal point position of the luminous flux exiting from the point 
5 of object height ri . 

An optimizing function (merit function) is defined as follows . 



(Optimizing Calculation) 

Optimizing calculation is performed by use of optical design 



QLO software. That is, the optimizing function F is minimized with 
U\ h 4/ hs and h 8 as variables. 



An example of the measurement result will be described. 
Specifications of the lens to be inspected are as follows. 



F = Azi 2 + Az 2 2 + Az 3 2 + 



[Example] 




Lens Outer Diameter: 0.25 mm<£ 



00 



Effective Radius r 0 = 0.125 mm 



Length Z = 1.001 mm (about P/2) 



Center Refractive Index n 0 



= 1.682 



Measurement was made on five lens samples to be inspected 



20 



(of the same lot number) . 



Measurement conditions are as follows. 



Parallel Striped Pattern: 100 line-pairs per mm 



Measurement Wavelength: A = 



654 nm 



Microscope Objective Lens: Plan 40 magnifying power (NA 



25 



= 0.65) 
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Microscope Eyepiece: 10 magnifying power 
Light Source: NA ^ 0.1 

When the average value of Af 0 was regarded as being in the 
paraxial focal point position, the value of the secondary index 
distribution coefficient _g was as follows. 

g = 3.113 /mm 

Fig. 5 shows curvature of field of five lens samples to 
be inspected (the relation between the distance from the optical 
axis of each lens sample and the meridional focal point position) . 
Measured data are plotted on the graph shown in Fig. 5. Variation 
in lens samples of the same lot number was so small that reproducible 
data were obtained. Fitting was performed by use of optical 
design software "Oslo Six" made by Sinclair Optics, Inc. in the 
U.S. on the basis the graph of curvature of field to thereby 
obtain the following values. 

h 4 = +1.639 

h 6 = -3.20 

h 8 = +2.57 

Curvature of field calculated on the basis of these index 
distribution coefficients is expressed as a solid line in Fig. 
5. As is obvious from Fig. 5, the calculated values coincide 
with the measured values well. It is confirmed from this fact 
that the higher-order index distribution coefficients h 4/ h 6 and 
h* can be obtained accurately and that the gradient index 
distribution of the rod lens can be measured. 



As described above, the invention provides a method in which 
higher-order index distribution coefficients are calculated back 
by a fitting process on the basis of measurement of curvature 
of field. The curvature of field can be measured accurately 
even in the case where the diameter of the rod lens is small. 
Accordingly, even in the case where the invention is used in 
an optical system requiring performance in diffraction limited, 
the gradient index distribution can be evaluated with sufficient 
accuracy . 

In addition, the method according to the invention need 
not use any laser light source. -There is also an advantage in 
that the wavelength for measurement can be selected relatively 
freely. 



